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Abstract. Magnetic null points can be located numerically in a potential 
field extrapolation or their average density can be estimated from the Fourier 
spectrum of a magnetogram. We use both methods to compute the null point 
density from a quiet Sun magnetogram made with Hinode's NFI and from mag- 
netograms from SOHO's MDI in both its high-resolution and low-resolution 
modes. All estimates of the super-chromospheric column density (z > 1.5 Mm) 
agree with one another and with the previous measurements: 3 x 10 null points 
per square Mm of solar surface. 



1. Null Point Density 

Line-of sight magnetograms, such as Figure[U show the quiet sun to be composed 
of a disordered mixture of positive and negative polarity. The coronal magnetic 
field above is expected to be complex and vanish at a numb er of points called 
null points (|Schriiver fc Titld l2002l:lLongcope & ParnelfeOQgT). These null points 



might serve as a site of heating (Has samlll992l: ICraig fc McClvmond Il993h or of 



coronal jets (jShibata et al.lll992l ; ICirtain et al.ll200 



ng 



The high resolution magnetogram from Hinode's NFI (Figure with 
0.16" x 0.16" pixels) reveals fine spatial structure, however, within a relatively 
small field of view: 136" x 96". The MDI full-disk magnetogram (Figure [lb) 
has lower resolution but from its larger field of view it is possible to extract the 
same region (labeled A) or a larger 300" x 300" square centered at disk center 
(B), similar to those used bv lLongcope Parnelll ( 20081 ). 



Potential fields are extrapolated from lower boundaries composed of the 
NFI magnetogram (Figure [T^,) and sub-field A of the MDI magnetogram. This 
region is near enough to disk center that no compensation is made for surface 
curvature — the line-of-sight field is used as the vertical field, B z (x,y,0). The 
ex trapolation yields field c omponents on Cartesian grid points. The algorithm 
of Havnes Parnell ( 20071 ) identifies each point at which tri-linear interpolation 



of the gridded field vectors vanish: magnetic null points. 
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Figure 1. Magnetograms of a section of quiet Sun from 19 Sept. 2008, 
(scaled between ±100 G). (a) From the MDI NFI instrument, (b) From 
SOHO/MDI in full-disk mode. Outlines show the field of view of (a) and of 
the 300" x 300" images used below. 



The heights of null-points directly found this way can be summarized as a 
density, pn(z), shown as the solid curves in Figure [2^. The null column, N z (z), 
shown in Figure [2)3 is the number of null points, per surface area, above a given 
height. For the surface area of the entire magnetogram (the right axis) the 
number of nulls above z ~ 8 Mm falls to only a handful, and the cumulative 
density suffers from statistical errors. 

Comparison reveals that additional spatial resolution yields more null points 
in the potential extrapolation from an ideal photospheric surface (z = 0), but the 
additional null points are, as expected, at very low altitude. Hence, the number 
of coronal null points {z > 1.5 Mm) is much less dependent on resolution. 



2. The Spectral Estimate 

An alternative to direct identification of null points is an estimate b ased on the 
isotropized Fourier spectrum of the magnetogram, S(k), developed bv lLongcope et al 



(2003). Figure [3^ shows the spectra from various magnetograms, computed out 



to k c = 2y/it/ 'Ax in order that the the total Fourier-space area, 7rk 2 , matches 
that occupied by the actual magnetogram with Ax x Ax pixels. 

The spectral estimate of null density depends on the variance, cr 2 (z), of the 
vertical component of the extrapolated field as well as on an inverse characteristic 
length scale, q(z), defined by a related integral 

roo 27T r°° 

a 2 z (z)=2ir S(k)e- 2kz kdk , q 2 ( z ) = —- S(k)e' 2kz k 3 dk . (1) 
JO °z\ z ) JO 



Longcope et al.l (|2003l ) showed that under the assumption that B z (x, y, 0) was a 
homogeneous, Gaussian random field, the density of magnetic null points is 

p N (z) = G(B z /a z )q\z) , (2) 



where G(£) is a dimensionless function ( Longcope et al. 120031 ). The function 



depends on the mean vertical field strength B z and also weakly on another 
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Figure 2. The density of null points above the Hinode magnetogram (red) 
and sub-field A of the MDI full disk magnetogram (blue), (a) The density, 
Pn{z), of null points found on a gridded extrapolation (solid) and by spectral 
estimate (dashed), (b) The column density, N n (z) found by integrating pn(z) 
downward. Left axis shows the value in Mm~ 2 , the right axis gives the total 
number of null points found on the grid. 



moment of the spectrum. The spectral estimate, Equation ([2]), is plotted in 
Figure [2^,, as a dashed curve for each of the two magnetograms. Integrating 
Pn{z) downward from infinity gives the column densities, N n (z), plotted as 
dashed curves in Figure [5)3. In the case of NFI it appears that the spectral 
estimate is reasonably accurate over heights with more than a handful of null 
points above. 

The spectral estimate depends on the isotropized spectra primarily through 
integrals in Equation (pQ). Different patches of the quiet Sun, such as regions 
A and B of Figure [U have very similar spectra, as shown in Figure Larger 
fields, such as B, sample S(k) more densely and thus yield more accurate spectra. 
Higher resolution, on the other hand allows spectra to extend to much higher 
wave numbers, as does the NFI curve in Figure^. The spectra from instruments 
of all resolution and fields-of-view agree relatively well over the lowest wave 
numbers, k < 1 rad/Mm. 

The moderate wave numbers of the MDI spectra are affected by white 
noise and the modulation transfer functio n of the imaging sy stem, including 
an intentional defocusing in full-disk mo de ( Scherrer et al.lll995l ). Correcting for 
these effects ( Longcope fc Parnelll I2008T ) yields spectra even closer to those of 
NFI. The NFI magnetogram is corrected by deconvolution with a PSF calculated 
from the lunar transit of July 2008; this has a small but noticeable effect in the 
data. 

Using each of the spectra from Figure [3^t in the spectral estimate leads 
to column densities shown in Figure [3)3. The fact that all of these agree over 
z > 700 km follows from the exponential factor in Equations ([T]). Due to this 
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Figure 3. (a) The isotropized spectra, S(k), for different magnetograms, 
plotted in different shades, offset vertically for clarity. Dashed lines are S(k) ~ 
e -2kd ^ wjth d — 12 Mm, for comparison, (b) Spectral estimates of the null 
column densities using the same shades. Triangles mark the column above 
z = 1.5 Mm. The dashed line is from the NFI spectrum, but with B z = 0. 



factor the spectral estimate at a height z will depend on the spectral range 
k < 1/z, where all spectra are similar. The fine-scale structure contributing 
to the high wavenumbers in both the Hinode NFI and the high-resolution MDI 
spectra create a layer of null points close to the photosphere (z < 500 km). 
Departures over z > 10 Mm are due to differences in k < 0.1 rad/Mm, where 
the different fields of view lead to different spectral forms. 

Evaluating the columns just above the chromosphere, at z = 1.5 Mm, gives 
JV n (1.5Mm) = 3.9 x 10~ 3 Mm" 2 for NFI and iV n (1.5Mm) = 3.1 x 10~ 3 Mm" 2 for 
MDI in high resolution (triangles in the inset). The latter is almost exactly the 
value found in a sample of 562 MDI quiet Sun magnetograms over 1996-1998 



and 2006-2007 (jLongcope fc ParnelJ |2Q0J ). Both the NFI and MDI full-disk 



estimates from field A are greater over the range 1 < z < 2 Mm (expanded in the 
inset) primarily due to the mean field over the small field of view, A: B z ~ —2.5 
G. The larger field of view (B) is much more flux-balanced {B z ~ —0.28 G). 
Using the NFI spectrum with B z = in Equation §2§ gives the estimates from 
the dashed line, falling much closer to those from field B. 
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